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A One-Dimensional  Fourier  Analogue  Computer 

Leonid  V.  Azaroff  and  M.J.  Buerger 
# 

Crystallography  Laboratory!  Massachusetts  Institute  of  Technology 


Abstract 

A computer  has  been  constructed  to  sum  Fourier 

m 

series  having  up  to  30  terms.  Although  this  is  a one- 
dixaensional  computer  it  can  be  used  for  double  and  triple 
summations  by  using  standard  trigonometric  expansions. 

4 

Secondarily,  it  can  be  used  for  computing  trial  structure 
factors • 

O.D.F.A.C.  sums  ZSa  qqq  2fftix  electrically.  The 

n — 

trigonometric  function  is  produced  by  a variable-angle 
transformer  known  as  a resolver.  Each  amplitude  is  set 
by  a variac  which  regulates  the  input  to  a particular  re- 
solver. The  frequencies  2 frhx  for  31' values  of  n are 
arranged  by  gearing  the  rotors  of  the  resolvers  in  ratios 

4 4 0 

0,  1,  2,  . . . 30.  The  resultant  individual  currents  are 
added  in  parallel  and  the  value  at  point  x (in  intervals  of 

4 4 

JL,  -1—,  or  1-  of  a cell  edge)  is  read  on  a voltmeter  and 
60  120  240 

the  phase  is  read  on  an  oscilloscope. 

The  relative  speed  of  a computation  is  five  to  ten 
times  faster  than  the  standard  strip  methods.  The  average 
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error  in  a computation  compares  favorably  with  the  rounding- 
off  error  in  conventional  tv/o-place  strips. 


Introduction 

Earlier^ computers . Fourier  series  are  important 
in  many  branches  of  science.  One-dimensional,  two-dimen- 
sional, and  three-dimensional  Fourier  series  are  especially 
important  in  x-ray  crystallography.  Since  the  computation 
of  these  functions  is  tedious,  a number  of  devices  have 
been  developed  to  perform  the  computation,  only  a few  of 

which  have  come  into  common  use.  Among  these  are  the 

l 2' 

electrical  digital  devioes  of  Beevers  * , the  electrical 

• » • 

analogue  machines  of  Hagg  and  Laurent^,  Ramsay,  et  al*f, 
the  mechanical  analogue  devices  of  MoLachlan  and  Champaygne^ 

6*7  g 

Rose  , Vand'  and  Beevers  and  Robertson  . All  these, devioes 

sum  one-dimensional  series.  Robert son^,  PepinskyJ‘°>11, 

and  McLaohlan  et  alXfc  have  devised  electrical  analogue 

maohines  for  two-dimensional  Fourier  summations, 

Basis^for^design.  When  the  number  of  Fourier 

syntheses  to  be  computed  in  the  Crystallographio  Laboratory 

of  became  large  enough  to  warrant  using  a special 

computing  devioe,  several  of  these  maohines  were  closely 

• 

investigated.  Sane  were  found  to  have  obvious  defects, 
suoh  as  oontaot  trouble  when  multiple  telephone  switohes 
were  used.  Most  of  them  were  found  to  handle  too  limited 
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a number  of  Fourier  terms..  Guided  by  this  survey,  a 

decision  was  reach.3i  to  build  a cne-dimenslonal  eleotri- 

oal  analogue  compute'*'  which  would  sum  a comparatively 

large  number  of  Fourier  terms.  This  was  specifically  set 
0 

at  30,  since  this  is  about  as  high  as  ever  required  in 
the  analysis  of  ordinary  non-protein  crystal  structures. 
McLachlan^  computer12  served  as  a guide  and  as  a point  of 
departure.  His  machine  compounds  phases  electrically  by 
using  selsyns.  It  performs  a two-dimensional  synthesis 
but  is  limited  to  8x8  terms, 

0 0 

The  transformation  of  an  angle,  y , into  a 
trigonometric  function  oan  be  performed  electrically  In 
many  different  ways.  The  two  devices  considered  in  the 
early  stages  of  designing  ODFAC  were  sine  potentiometers 
and  resolvers,  A sine  potentiometer  consists  of  a con- 
tinuous resistor  winding  tapped  to  produce  a voltage  which 
varies  as  the  sine  of  the  angle  of  rotation  of  the  main 
shaft.  A i-osolver  is  a transformer  U3ing  as  its  primary 
the  rotor  member  and  as  its  secondary  two  separate  stator 
windings  placed  90  degrees  apart.  When  a voltage  is 
applied  to  the  primary  winding  the  voltages  of  the  two 
secondary  windings  vary  respectively  as  the  sine  and  co- 
sine of  the  angle  of  rotation*  Resolvers  were  chosen  for 
the  construction  of  ODFAC  bocauso  the  isolating  property 
of  a transformer  makes  the  adding  circuit  independent  of 
variations  in  the  input  stages. 
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Principle  of  operation.  Fig.  1 shows,  in  outline 

0 

form,  the  operation  of  QDFAC.  The  input  to  each  resolver, 

4 0 

R,  is  fed  from  a variao,  V.  The  variao  thus  controls  the 
amplitude  or  Fourier  coefficient,  A,  of  the  Fourier  compon- 
ent of  a particular  resolver.  The  two  outlet  leads  of  the 
resolver  then  deliver  voltages  proportional  to  A cos  <f>  and 
A sin  <p  . The  shafts  of  30  such  resolvers  are  geared  so 
that  the  shaft  displacement  of  a particular  resolver,  n, 

0 0 

is  an  integral  multiple,  n,  of  the  shaft  displacement  of  a 
fundamental  shaft.  The  outputs  of  the  individual  resolvers 
are  therefore  a set  of  voltages  proportional  to 

A0 

A^  cos  <f±  sin 

Ag  cos  if 2 A-2  sin  <f2 


A^qCos  A^QSin  fjQ 

If  the  outputs  are  appropriately  coupled,  the  machine  pro- 
duces voltages  proportional  to 

30  30 

V.  An  008  /n  “A  E An  sin  yn 

n-o  n-o 

which  can  be  read  on  a voltmeter,  VM,  In  orystallographio 
problems  - 2m  where  n is  the  number  of  the  harmonio 
and  x is  the  sampling  interval  expressed  as  a fraction  of 
one  complete  period. 
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Design  and  construction 

Mechanical  features.  The  product t on  of  thirty 
harmonics  in  the  Fourier  series  requires  that  one  complete 
rotation  of  the  rotor  in  the  first  resolver  corresponds 

to  two  oomplete  rotations  of  the  rotor  in  the  second  re- 

* • 

solver,  etc.,  up  to  thirty  complete  rotations  of  the  rotor 
in  the  thirtieth  resolver.  This  is  accomplished  in  ODFAC 
by  means  of  a gear  train  represented  by  Fig.  2, 

The  main  shaft  leading  from  the  motor  is  geared 
by  worms  and  worm  gears  to  three  horizontal  shafts  whose 
relative  angles  of  rotation  are  given  by  the  ratios  of  the 
worms  and  gears.  Each  horizontal  shaft  carries  a series 
of  spur  gears  which  engage  with  the  spur  gears  mounted  on 
the  rotors  of  the  individual  resolvers.  The  angle  of  ro- 
tation of  each  rotor  shaft  is,  therefore,  a function  of 
the  ratio  of  its  spur  gears  modified  by  the  angular  rota- 
tion of  the  horizontal  shafts.  The  actual  gear  train  is 
shown  in  Fig.  3.  This  particular  gear  train  has  several 
advantages.  One  is  that  most  gears  used  are  commercially 
available,  stock  sizes.  Only  10  of  the  82  gears  are  non- 
standard size  and  had  to  be  obtained  by  special  order. 
Another  advantage  i3  that  high  gear  ratios  are  avoided, 
the  largest  step-down  ratio  being  5:4»  This  decreases 

inaoouraoies  in  the  positioning  of  the  resolver  shafts • 

0 

Because  of  low  friotional  losses,  all  30  resolvers  can  be 

0 

driven  at  the  desired  speed  by  a small,  hip.  motor. 
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Electrical  circuit*  The  wiring  diagram  of  the 
eleotrioal  components  of  ODFAC  is  shown  in  Fig,  4*  The 

4 4 

input  voltage  is  fed  from  the  variacs,  V,  through  a DPDT 

• • 4 

phase-selector  switch,  S,  to  the  rotors  of  the  resolvers, 

4 * 

R,  which  are  linked  by  the  gear  train  illustrated  in  Fig, 

3,  Each  output  line  of  the  stator  winding  of  the  resolvers 

4 » 4 

contains  a series  resistance,  r,  of  500,000  ohms  whose 
purpose  is  to  make  any  stray  losses  in  the  secondary  cir- 
cuit negligible*  All  the  cosine  (and  sine)  windings  of 
even  harmonics  are  connected  in  parallel.  Similarly 
the  cosine  (and  sJjie)  windings  of  the  odd  harmonics  are 
connected  in  parallel.  These  lines  lead  to  a gang  switch 
which  permits  various  combinations  of  the  lines  to  be  made. 
The  use  of  a parallel,  rather  than  a series  adding 
circuit  is  preferred  for  the  following  reasons:  In  a series 

circuit  all  the  stater  windings  would  be  directly  connected 
causing  a cumulation  of  errors  due  to  mutual  inductance 
between  the  individual  primary  and  secondary  windings  and 
the  secondary  windings  of  adjacent  resolvers.  The  addition 
of  voltages  in  a series  circuit  would  also  have  the  dis- 
advantage of  building  up  very  high  voltages  and  currents 
with  possible  damage  to  the  windings.  All  of  these  dis- 
advantages are  overcome  by  a parallel  circuit  in  which 
the  isolation  of  the  individual  resolvers  limits  tho  ourrent 
in  each  secondary  winding  to  that  induced  by  the  primary. 
This  isolation  permits  tho  insertion  of  a large  resistance 
to  minimize  losses  in  the  line.  The  maximum  output  voltage 
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obtained  from  suoh  a circuit  oan  never  oxoeed  the  input 
line  voltage  beoause  the  individual  voltages  are  averaged 
over  all  31  branohes. 

The  output  voltage  lines  pass  through  a 4-gang 

I 

SWitoh  which  selects  combinations  of  odd  and  even  harmonios  ! 

i 

of  the  sine  or  oosino  lines  for  transmission  to  the  volt- 
meter and  oscilloscope,  The  separation  of  the  odd  and  even 
harmonics  into  separate  lines  permits  the  utilization  of  the 
symmetry  inherent  in  sine  and  cosine  functions*  This  is 

standard  practice  in  orystallographic  applications1^  The  1 

? 

meter  has  a 50/* a movement  and  three  ranges  which  oorre- 

spond  to  the  voltages  expected  from  the  number  of  variaos 

* 

supplying  input  voltages.  A fourth  range,  actuated  by  a 

push-button,  permits  a clear  reading  of  low  voltages. 

The  oathode-ray  tube  is  part  of  a standard 

osoillosoope  oiroult  not  having  a sweep  generator.  The 

output  voltage  is  placed  aoross  two  deflecting  plates  and 

a reference  voltage  aoross  the  other  two.  Since  both  sets 

of  plates  have  the  same  frequency  applied  to  them,  the 

resulting  Lissajous  figure  on  the  tuba  face  is  a straight 

line  whose  angular  inclination  depends  on  the  magnitude  and 

relative  phase  of  the  output  voltage « 

Constructional  details.  The  rear  view  of  the 

assembled  machine  is  shown  in  Fig,  5*  Eh©  resolvers  are 
* # 

held  in  place  by  aluminum  clamps  (see  also  Fig.  3)  which, 
in  turn,  are  fastened  to  two  aluminum  angles  forming  a 


I 

> ... 


e9 


track-like  shelf,  Each  shelf  holds  five  resolvers  and 
supports  its  horizontal  drive  shaft  mounted  in  ball- 
bearing supports.  The  clamps  have  built-in  means  for 
adjusting  the  angular  position  of  the  resolvers.  The 

horizontal  drive  shafts  are  driven  by  the  main  drive 

* * 

shaft  which  runs  vertically,  in  the  center  of  Fig,  5, 
from  the  motor  (hidden  by  the  bottom  shelf)  to  the  con- 
trol mechanism. 

The  control  mechanism  is  connected  directly  to 
the  drive  shaft  by  a set  of  ohangG  gear3.  It  consists  of 
a shaft  which  carries  a change  gear  and  a cam*  The  cam 
actuates  a micro-switch  which  controls  the  number  of 
revolutions  made  by  the  motor  for  the  sampling  interval. 

The  motor  is  a dynamic -braking  motor  which  can  be  stopped 
instantaneously  by  applying  a reverse  field  to  its  rotor* 

A horizontal  shaft,  permanently  geared  to  the  main  drive 

4 

shaft,  rotates  a dial  on  the  front  panel  indicating  the 
angular  position  of  the  main  shaft.  The  change  gears  in 
the  control  mechanism  permit  the  selection  of  intervals 

4 4 

of  JL,  JL,  or  -i-  of  one  complete  period, 

60*  120  240 

The  motor-generator  set  at  the  bottom  of  Fig,  5 
provides  the  input  voltage  (120  v,,  400  ~)  to  the  varlaos 
only.  The  driving  meohenism  and  service  components  operate 
on  the  regular  line  voltage.  The  purpose  of  the  motor- 
generator  sot  is  to  provide  a oonstant  current  souroe 

4 

with  an  undlstorted  wave  shape,  making  the  adding  oirouita 
in  0DFAC  independent  of  line  fluctuations* 
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Performance 

Operation#  A front  view  o.?  CPFAC  in  operation 

ie  shown  in  Fig.  6#  The  variaos  ere  numbered  and  arranged 

on  removable  panels  in  sets  of  five#  The  amplitudes  are 

set  by  means  of  a frict ion-drive  dial  permitting  an 
# 

acourate,  rapid  setting  from  0-100c  The  phases  are  set 
by  means  of  a plus-minus  toggle  switch  placed  direotly 
above  each  dial.  Eaoh  banlc  of  five  resolvers  can  be 
switohed  in  or  out  of  the  input  oirouit  by  means  of  an 
additional  toggle  awitoh  on  the  left  of  eaoh  panel*  The 
numbering  and  looation  of  the  dials  is  such  that  the 
operator  sits  faoing  the  first  sixteen  dials,  the  next 
fifteen  dials  being  within  arm's-length  to  hie  right* 

The  dials  need  be  set  only  onoe  for  eaoh  one-dimenalonal 
summation* 

A panel  to  the  left  of  the  operator  contains 
the  meter  on  whioh  the  value  of  the  Fourier  aeries  is  read 
the  CRO  tube  on  whioh  the  positive  and  negative  quadrants 
are  marked,  the  sine-oosine  selector  switch,  and  a push- 
button that  advances  the  computer  mechanism  to  the  next 
sample  setting.  A pilot  light  indicates  that  the  motor 
has  advanced  the  computer  to  the  n*xt  sample  setting  and 
the  meter  may  be  ror.d*  Onoe  the  amplitudes  (variao  dials) 
and  the  phases  (pliu!~winus  switches)  have  been  set,  and 

the  type  series  (sine  or  oosine)  has  been  seleoted,  all 

• # 

the  operator  does  is  read  the  meter,  note  the  phase,  and 
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push  the  button  to  advance  to  the  next  reading.  A drop- 
leaf  table  is  attaohed  to  the  front  of  0D!?AC  to  provide  a 
convenient  surfaoe  for  recording  readings. 

The  speed  of  a computation  on  ODFAC  depends 

* * 

primarily  on  the  type  of  series  deoired,  i.e,,  one- 

* * • 

dimensional,  two-dimensional,  eto*,  and  only  seoondarily 
on  the  number  of  terms  or  frequency  of  interval  of  sampling 
desired.  The  maohine  time  for  a oomplete  oyole  from  0 to 
2 ir  has  been  seleoted  to  take  7 minutes.  If  the  time  for 

setting  the  variaos  and  recording  the  values  obtained  is 

# 

added  to  this,  the  total  computing  time  amounts  to  approx- 
imately  10-15  minutes,  depending  on  the  sampling  interval 
selected.  The  relative  time  of  computing  an  average  two- 
dimensional  series  on  ODFAC  ranges  from,^  to  -L  of  the  time 
consumed  using  standard  strip  methods1*1, The  more  terms 

there  are  in  the  series  and  the  finer  the  interval  of 

« 

sampling  desired,  the  more  efficient  ODFAC  beoomes  when 
compared  with  strip  methods. 

Aoouraoy.  The  aoouraoy  of  the  components  used 
in  ODFAC  is  the  highest  attainable  at  a reasonable  oost. 

The  resolvers  are  aoourate  to  within  lj^  meohanloal  degrees. 
The  variaos  and  resistors  are  aoourate  to  within  1%  of 
maximum  ratings.  The  meter  is  aoourate  to  better  than 
1%  of  full  soale  defleotlon.  The  baoklash  in  the  gears 
is  almost  non-existent  and  the  angular  aoouraoy  of  the 
gear  settings  is  within  a fraotion  of  one  degree. 

Table  1 shows  an  aotual  comparison  of  the  meter 
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readings  of  ODFAC  with  those  oomputod  using  three-place 
trlgonooetrio  tables.  The  amplitudes  for  this  one- 
dimensional  series  were  supplied  by  Professor  M. J.  Buerger 
from  a Barker  line  synthesis  for  realgar.  In  the  same 
table  are  listed  corresponding  values  as  oomputed  with 
the  aid  of  Patterson-Tunell  strips1**'  and  Beevers-Lipson 
strips1 If  the  deviations  from  the  true  values  ard 
examined  it  is  evident  that  the  errors  in  the  values 

• 0 9 0 

given  by  ODFAC  are,  on  the  average,  as  small. as,  if  not 

• 9 

smaller  than,  those  due  to  the  rounding-off  errors 

inherent  in  the  strip  methods « 

* 

Finally,  by  making  all  operations  other  than 

0 

recording  the  numerical  values  automatic,  OQFAC  eliminates 
the  ever-present  souroe  of  error  - the  human  error. 

OohoiUsions  * Tke  thief  advantages  of  ODFAC 
are  its  ability  to  handle  up  to  31  coefficients  in  a 
Fourier  series  and  to  perform  the  summation  rapidly  with 
an  aoouraoy  adequate  for  most  purposes.  The  simplicity 
of  the  design  of  the  maohine  virtually  eliminates  almost 

all  possibilities  of  eleotrioal  or  meohanlcal  failure. 

* 

In  the  event  suoh  failure  should  ooour,  all  components 
are  readily  aooesslble  and  removable  for  repair. 

The  slowest  part  of  the  operation  of  ODFAC 
lies  in  recording  the  results.  It  is  proposed  to  add  an 
automatic  reoorder  for  this  purpose* 


Extension  to  two  dimensions 
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Use^of ^two^unii s , One  CJPAC  unit  performs  a 
two-dimensional  synthesis  by  meacs  of  a Esevers-Lipson 
expansion*  A combination  of  two  such  units  performs 
this  summation  more  rapidly.  If  these  units  are  equipped 
with  ohart-type  recorders,  the  most  complicated  two- 
dimensional  series  having  ?0  x 20  Courier  coefficients 
can  be  performed,  in  an  afternoon.  The  utilization  of 
two  such  units  follov/s  the  plan  Indicated  below. 

The  series  to  be  summed  has  a general  expression: 

fi(xf)  • I 5?*hk  609  2?f(hx+ky) 

[ h k 

■ ri[^Fhk ooa  27}*y]  ooa  2frbx  • ^f£yhk  sln  2wi&] 9in  2*hx 

The  summations  over  k are  first  performed  at  the  same 
time  by  the  two  units  working  independently.  The  results 
of  these  first  summations  are  then  used  as  ooeffiolents 
for  the  seoond  summation  over  h.  In  the  second  operation 

0 

the  two  units  are  oonneoted  in  parallel  t.o  only  one  reoorder, 
whioh  thus  records  the  results  of  the  final  two-dimensional 
summation. 

^o-dimenflional^ocmjuter  • A more  elaborate  ex- 
tension to  two-dimensions  can  be  built  by  combining  two 
resolvers  for  every  term  in  the  series.  The  oonstruotlon 
entails  building  k ODTAC  units  and  one  "master"  ODFAC 
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having  h resolvers  (without  variaos)  for  each  harmonio  k. 

Suoh  a plan  is  illustrated  in  Fig*  7 for  k » h - 4. 

The  ODFAO  units,  whose  resolvers  have  angular 

0 

Speed  ratios  proportional  to  their  particular  h,  are 
Shown  in  solid  lines.  The  reotangles  in  the  dotted 
lines  represent  the  resolvers  of  the  ’’master"  unit. 

All  Of  these  latter  resolvers  have  the  same  angular 
settings  proportional  to  k. 

The  coupling  between  all  units  is  electrical 
only.  The  resolvers  of  one  ODFAC  unit  produce  voltages 
proportional  to: 


V* 008  2*v 

Fhjli}  008 

«?8 


which  are  fed  to  the  resolverB,  R^,  of  the  master  unit. 
The  outputs  of  these  resolvers  are  thus  proportional  to: 

^*1*1  oos  2*&ix)  008  2rkiy 
^yh2k1  008  008  2^7  etc. 

Appropriate  eleotrioal  coupling  of  the  resolvers  in  the 
master  unit  then  produces  a voltage  proportional  to 

£ £ *hk  008  2rr(hx  ♦ ky) 


The  authors  wish  to  express  their  gratitude 
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of  the  computer. 
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Ono-Dimanr: icuai  Jourie?  Synthesis 
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Fig.  3 Gear  Train 


Control  Mechanism 


Fig.  5 ODFAC,  Rear  View 


Fig.  6 ODFAC,  Front  View 
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